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ABSTRACT 


Cyclotrimethlenetrinitramine (RDX or Cyclonite) exists in two poly- 
morphic species, a and 8. ‘The a-form is stable under normal conditions, 
whereas 8-RDX is formed by recrystallization from high boiling solvents. 
It immediately transfers to the a-form in the presence of excess a-RDX, 
and itis therefore assumed that"RDX as it exists in its normal state 
consists solely of the a polymorph. 

The thermal decomposition of RDX has been comprehensively studied 
at Picatinny Arsenal by means of quantitative measurements of the de- 
composition products. The thermal stability of explosive compositions 
has also been investigated by measuring the rate of gas evolution during 
decomposition. 

It was the purpose of this work to obtain thermal and kinetic data 
for RDX using Differential Thermal Analysis (DTA) experimental methods. 
These results and techniques were compared with those obtained using 
the other methods and with theoretical values. The favorable results 
lend credibility to the DTA method and suggest that further development 
in the technique will prove useful in providing a fast and reasonably 
accurate method of characterizing the thermal and kinetic properties of 


new explosives and propellants. 
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I. INTRODUCTION 


Explosives undergo thermal deconposition at temperatures far below 
those at which spontaneous explosion may oceur. The decomposition re- 
action is important in determining the stability of the explosive. А 
considerable amount of heat energy is liberated in the decomposition 
process. This highly exothermic reaction aecelerates to a high veloc- 
ity culminating in deflagration or detonation when the rate of heat 
production is greater than the rate of heat loss to the surroundings 
acswietion amdseonveetjon. Explosives in their normal state are 
metastable, the free energy of the explosive being considerably higher 
than its decomposition products. For an explosive to remain stable at 
ordinary temperatures, its decomposition reaction must depend on rel- 
ОШ Iy high activation energies. 

Kinetic studies are therefore important in the study of explosive 
reaction mechanisms. The Differential Thermal Analysis technique is 
one of several that is used to study reaction kinetics. When a reac- 
tion occurs in DTA, the change in heat content and in the thermal prop- 
erties of the sample is indicated by a deflection in the thermogram. 
For reactions that possess activation energies and thereby proceed at 
a rate varying with temperature, the position of the peak varies with 
the heating rate. The activation energy and Arrhenius frequency factor 
can be determined from this experimental information. 

Cyclotrimethylenetrinitramine (RDX or Cyclonite) is an important ex- 
plosive used in many military composite explosives. RDX melts at 203°C, 


and the decomposition of liquid cyclonite follows the unimolecular rate 





equation over the range of temperatures 213-299°C. [Robertson 1949]. 
The Arrhenius frequency factor is large, and the activation ‘energy is 
of the order of 50 kcal/mole. The high value of the Arrhenius frequency 
factor is a possible indication of self heating. 

The purpose of this work was to investigate the thermal decomposi- 
tion of RDX using DTA experimental techniques. In the decomposition 
of 1,3,5,7-Tetranitro-1,3,5,7-Tetrazacyclooctane (HMX), three separate 
peaks were observed in the DIA thermogram, and three separate activa- 
tion energies obtained from the DTA data [Нопдее 1971). It was of 
interest to determine whether these phenomena were evident in the de- 
composition of RDX, to detect any evidence of self-heating, and to 


compare the kinetic data obtained with that found in the literature. 





ІІ. THERMAL DECOMPOSITION OF EXPLOSIVES 


The investigation of the thermal behavior of explosives oes both 
the study of chemical kinetics and the theory of transport properties. 
The chemical kinetics of explosive reactions are not fully understood due 
to the high reaction rates and temperatures involved and the complexity 
ocwbhe"peodnets of "therreaction. "Therdecomposttion-of RDX may be written: 


O ల [+ ος 
3666 2 2 2 2 


This reaction may be accurate for detonation but is unlikely to occur 
thermally due to the large number of bonds which must be broken in the 
parent molecule. Attempts to study these reactions require that mea- 
surements of the decomposition rates be made at relatively low 
temperatures. 

Thermal decomposition is essentially a heat-balance problem. The 
heat loss to the surroundings and the accumulation of heat in the ex- 
plosive equals the chemical energy generated by the decomposition of 
сре explosive. 

KV?T + pcd = da/dt 

This can be written simply as F + G =H. The chemical energy generated 
by the decomposition, H, increases exponentially with temperature. The 
heat loss, F, consists of thermal conduction, which is proportional to 
the thermal gradient, and radiation, which follows the Stefan-Boltzmann 
T* law. The heat loss therefore increases at a slower rate, especlally 
since the heats of activation of explosives are of the order of 50 
kcal/mole. The accumulation of heat in the explosive, G, becomes more 
important as the temperature increases, causing acceleration of the 


decomposition. 
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Some explosives liberate energy in their interior due to a slow chem- 
ical reaction that is often a decomposition reaction. This is known as 
self-heating and causes the temperature of the material to increase, re- 
sulting іп an increase in the rate of reaction. Dynamic equilibrium is 
reached when the heat is removed as fast as it is generated. If heat is 
generated faster than it ean be removed, the temperature increases more 
rapidly and the reaction accelerates. The maximum temperature of the 
surroundings for which a steady state is possible is called the critical 
temperature for the explosive. Above this temperature, the material will 
self-heat until ignition or explosion occurs [Longwell 1961]. 

Decomposition reactions are sometimes autocatalytic or catalyzed by 
small amounts of impurities. The decomposition products or impurities 
ete los ive nay act as a catalyst in that the reaction rate is in- 
creased as more products are produced. The decomposition rates of these 
materials then become a function of the duration or extent of the 


decomposition. 


ІШЕ 





III. REACTION KINETICS IN DTA 


The material of this section is a summary of Kissinger's work on 
the theory of differential thermal analysis [Ref. 2]. DTA techniques 
produce a deflection or peak (AT) in the thermogram due to a change in 
heat content and thermal properties of the active sample during a re- 
action. This peak temperature is affected by experimental technique, 
and is often higher than known transition or decomposition temperatures. 

The temperature of maximum deflection is defined by kinetic con- 
stants and the heating rate. If a reaction posses an activation energy, 
it will proceed at a rate varying with temperature and the position of 
the peak will vary with the heating rate. The temperature of maximum 
deflection in DIA is also the temperature at which the reaction rate is 
a maximum. The temperature distribution of the sample explosive obeys 
the general heat flow equation; 


ot pc pe dt 


Diez cher rate of heat generated due to chemical reaction per unit vol- 
ume of the sample is equal to zero in the inert reference. The temper- 
ature distribution in the reference material is then given by; 


2 


с pc 
Assuming the sample to be a cylinder of radius a and of infinite 
length, integration of equation 2 yields; 
| Pe Ty et — Spent (3) 
where ф - dT/dt and the outside temperature is given by T = To + $t. 


The rate of heat generated due to chemical reaction is a function of 
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temperature in the active sample. Equation l is therefore a non-linear 
partial differential equation. Assuming that the temperature of the out- 
side of the tube rises at a linear rate, the solution for the temperature 
at the center of the sample will be of the form; 

Ts = Th + ét - f(dq/dt) (4) 
where f(dg/dt) is a function of the reaction rate and includes secondary 
effects of the reaction such as changes in volume, density and thermal 
properties. 

The differential temperature is the difference between temperatures 


ne" cembers of the two samples; 





" (нее; 
к (геГерепсе (5) 


= г(28) арте 
and is a maximum at the temperature of maximum deflection where de/dt 
is zero. 

de/dt = f'(dg/dt) d?q/at? (6) 
When the derivative of the rate of heat absorption is zero it can be 
seen from Equation 6 that de/dt is also zero. Since the rate of heat 
absorption is proportional to the rate of reaction, Equation 6 shows 
that the peak differential deflection occurs when the reaction rate is 
a maximum. It should be noted that the heating rate must be constant 
for this to hold and that Equation 6 is valid for a sample of any shape. 
The rate law for the thermal decomposition of explosives can be 


described by the equation; 
(3x/3t)p 7» kp (1-x)H (T) 


where km = A exp(-Ea/RT) and temperature is a function of time. 


dx . (ax) , fax|, ат 8 
` at (æt) + (57) dt p 
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The rate of change of x with respect to temperature, with the time 
coordinates fixed, is equal to zero because the number and positions 
of the particles of the reactant are also fixed with time. The only 
effect of instantaneous change in temperature is in the velocities of 
thermal motion of the reactant particles. Equation 7 can now be 
written; 

dx/dt - A(1-x)"exp(-Eg/RT) (9) 
which holds for any temperature as long as the fraction of reactant 
decomposed and temperature are measured at the same instant. ЈЕ 15 
assumed that the empirical order of the reaction remains constant 
πα ле ereater part of the reaction. 

When the temperature rises during the reaction, the reaction rate 
will rise to a maximum “value, then return to zero as the reactant 15 
exhausted. This maximum rate occurs when d/dt(dx/dt) is zero. Pro- 
ο ο σας егатше rises at a constant rate, the differentiation 


of Equation 9 yields; 


а (ах) _ ах | Еаф -l „-Ед/ВТ 
Ты n. o -x)n 
а (= at Ë: An(l-x) = e78 (10) 
Setting this equation equal to zero defines the temperature at which 


the maximum rate occurs. 


RT? 


= дп(1-х)071 е7Еа/ (11) 
IE. ss the затрје temperature at which the peak in the DTA thermogram 
occurs. The quantity of material left unreacted, (1-x),,, is not de- 
termined by the DTA method. 


To obtain the extent of reaction as a function of temperature, 


Equation 9 can be integrated. The resulting exponential integral does 
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not yield a simple expression. Previous work has resulted in a satis- 
factory approximation by successive integration by parts. Higher order 


terms were neglected, and the simplified expression becomes; 


п(1-х)й-—1 - ] + (n-1) = (ТО) 


The heating rate does not appear except as a result of T,,, which varies 
with heating rate. The product п(1-х)0-1 is independent of heating rate 


and nearly equal to unity. Equation ll now becomes; 


-Ea /RT 
9 = AR а/ πι (13) 
Еа 
Taking the natural logarithm of both sides and differentiating; 
d In($/T%) _ d In (AR/Ea eHa/RIm) 
a(1/Tm) d(1/Tm 
E, Ve*Ea/RIm (два (Lg, e ER 
ae 
2 -Е 
d In($/Tý) “Га (14) 


aooo) R 


By experimentally obtaining DTA thermograms at different heating rates, 
the activation energy for the decomposition reaction can be determined 


by Equation 14, regardless of reaction order. 
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IV. EXPERIMENTAL PROCEDURE 


A. APPARATUS 

The differential thermal analysis reaction was obtained by use of 
the DuPont 900 Differential Thermal Analyzer shown in Figure l. А ге- 
mote cell unit, shown in Figure 2, was used in this work due to the 
nature of the samples being studied. 

DTA uses a differential thermocouple arrangement consisting of two 


thermocouples wired in opposition as shown in Figure 3. 


SAMPLE 


REFERENCE 





FIGURE 3 DIFFERENTIAL THERMOCOUPLE 
Thermocouple A is placed in the sample to be analyzed. Thermocouple B 
is placed in an inert reference material. When the temperature of the 
sample equals the temperature of the reference, the thermocouples pro- 
duce identical voltage, and the net voltage output, AT, is zero. When 
a physical or chemical reaction occurs in the sample, a differential 
signal is received and recorded. Figure 4 shows a schematic of the DTA 
system. 

A silver heating block is mounted in the remote cell assembly. 
This block contains the sample, reference, and controlling chromel- 
alumel thermocouples as shown in Figure 5. The block is designed to 
provide for equal heating rates for both the reference and sample ma- 


terials. The temperature of the reference is less than the temperature 
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FIGURE 2: Remote Cell Assembly 
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FIGURE 4: DTA System Schematic 


(taken from ref. 10, pg. 3-2) 
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FIGURE 5: Heating Block and Thermocouple Arrangement 


(taken from ref. 10, pg. 3-4) 
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of the block by a constant amount that is directly proportional to the 
heating rate, weight and heat capacity of the reference material, and 

inversely proportional to the thermal conductivity of the reference ma- 
terial, plus a — exponential funetion which can be neglected for 


simplification. 
I To tp (We/k JR 


Assuming a thermally inert sample, 


Те=Тв-Ф (Ме/к) 


Т = temperature к = thermal conductivity 
þ = heating rate S = sample 

W = weight R = reference 

c = heat capacity B = block 

The differential temperature 15: 


АТ-Те-Ть-фр(МС/к)р-фе(МС/к)с 
౯-౫౯ block in the DuPont 900 is designed such that $s76p: 
AT= (We/k )R-(We/K )s 
If the weight, heat capacity and thermal conductivity of the sample and 
reference material are the same, AT is zero. The AT signal as a func- 
tion of sample temperature is sent to the X-Y recorder. Any change іп 
the properties of the sample results in a positive or negative value 


for AT as shown in Figure 6. 


POSITIVE SIGNAL: (We/k )p> (We/k )g 
e ,kz f (T) 
ZERO SIGNAL: (We/k )p=(We/k )g 


ENDO «АТ» EXO 


NEGATIVE SIGNAL: — (Wc/k)p<(We/k)g 


с,к#Р(Т) 
SAMPLE TEMP. > 


FIGURE 6 DTA THERMOGRAM 
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A thermally active material will experience a change in heat capac- 
ity and thermal conductivity over a temperature range. The AT trace 
will be a sloping straight line if the properties change linearly with 
temperature, КОЗ согуга ШІ спе change with temperature is quadratic or 
of higher order. Figure 7 exhibits this relationship, where the curved 


line would have a positive slope if (Ме/к )р> (Ме/к) с. 


um (We /k )y<(We/k)g 
| 1. linear f(T) 
(We/x )p<(We/k )g 

M quadratic f(T) 


ΕΝΡΟ «ΔΙ» EXO 


SAMPLE TEMP + 
FIGURE 7 DTA THERMOGRAM 
The DTA thermogram of a thermally active sample will exhibit endo- 
therms and exotherms where the temperatures of interest are defined in 


Figure 8. 


| 
| EXTRAPOLATED ONSET 
| 








ENDOTHERM «АТ»  EXOTHERM 


SAMPLE TEMPERATURE > 


FIGURE 8 TRANSITION TEMPERATURES 
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TRANSITION TEMPERATURE DEFINITIONS USED IN FIGURE 8. 


ONSET: The temperature at which the thermogram starts to depart from 
the Pase line: 


EXTRAPOLATED ONSET: Тһе temperature corresponding to the intersection 
of extrapolations of the base line amd the longest straight-line 
section of the low-temperature side of the peak. 

РЕАК: The temperature of reversal. 


RECOVERY: The temperature at which the thermogram returns to either 
the same or a different base line. 


EQUILIBRIUM: The region where the state or form of the sample before 
transition is in equilibrium with its state or form after 
transition. 

For the purposes of this work, a fourth well was drilled in the 
silver heating blocks to cross check the temperature scale on the DTA 
900 X-Y plotter. A single chromel-alumel thermocouple was set in inert 
glass beads and connected to a Varian recorder. This gave a millivolt 
versus time chart to check the heating block temperatures and heating 
rate. The Varian thermocouple was calibrated with a millivolt poten- 
tiometer and temperatures were verified through the use of a hot oil 
bath. The temperature scale of the DTA thermogram is calibrated to 
the chromel-alumel millivolt output to give direct temperature readings 


in degrees centigrade. 


В. SAMPLE PREPARATION 

The samples of RDX that were used in the experimental runs were 
prepared by drying a small quantity of pure RDX that was stored in 
water. The sample appeared to be pure RDX as its melting temperature 
of 204°C compared favorably with listed values. 

Two separate heating blocks were used in the experimental work. In 


-the variable heating rate runs, a silver block with wells to accept two 
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millimeter diameter capillary tubes was used. The RDX sample was 
tapped down into one of the four capillary tubes to a depth of four 
millimeters, which was иа to be 19.3 milligrams of RDX. The 
three remaining capillary tubes contained equal levels of inert finely 
ground glass beads .05 millimeters in diameter. These tubes contained 
the reference, controlling, and Varian recorder chromel-alumel thermo- 
couples. All four capillary tubes were then cut to 1-1/8 ineh. Care 
was taken when placing the = и to insure that they were 
pushed to the full depth of the material. The capillary tubes were 
then placed in their respective wells, and the thermocouple leads were 
folded under the rim of the remote cell unit to insure that they re- 
mained completely seated throughout the run. This was especially im 
portant in Che case of the RDX sample since Che reactant gases given 
off during the exothermie decomposition tended to push the thermo- 
couple out of the tube. The sample depth of four millimeters was 
critical for satisfactory results. A quantity larger than this would 
push the thermocouple out of the tube, giving erroneous readings. A 
smaller quantity would not provide complete coverage of the thermo- 
couple probe. The fineness or texture of the sample in this work was 
not considered since the reactions occurred above the melting tempera- 
ture, or in all cases for the variable heating rate runs, in the liquid 
phase. 

The second silver heating block was identical to the first with the 
exception of the sample and reference wells which were designed to hold 
four millimeter diameter test tubes. The RDX sample was measured in 


the same manner as before to a depth of four millimeters, which in 
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this case gave a sample of 43.32 milligrams of RDX. The second four 
millimeter test tube contained an equal volume of inert glass beads and 
served as the reference. Both four millimeter tubes were covered with 
small metal caps through which the thermocouples were inserted. This 
reduced losses to the atmosphere during reaction. The controlling and 
Varian thermocouples were inserted in the remaining two small capillary 
wells in the same manner as in the variable heating rate block. This 
second heating block was used in the isothermal portion of the 


experimental work. 


C. VARIABLE HEATING RATE EXPERIMENT 

The DIA 900 control panel contains dials to set starting tempera- 
tures and heating rates. It was determined that there was no apparent 
difference in the thermogram presentation for a sample that was heated 
from room temperature at an established rate through decomposition or 
for one which was heated at an accelerated rate (50°C/min.) to 200°C, 
then slowed to the established rate through the decomposition reaction. 
In order to expedite the runs, this latter procedure was adopted. 

The full range of the DTA 900 heating rate capability was utilized. 
In order to establish consistency in the procedure, an outline was set 
down and followed in each run. 

EXPERIMENTAL PROCEDURE FOR VARIABLE HEATING RATE EXPERIMENT 

(1) Turn DTA 900 power knob to STANDBY, allow five minute warm up. 


(2) Place inert reference and sample material in the capillary tubes, 
tap and measure to a depth of four millimeters. Cut to 1-1/8 inch. 


(3) Insert corresponding thermocouples, insuring that the thermocouple 


probe is completely inmersed in the material. A new thermocouple 
is used for the RDX sample in each run. | 
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(Ң) Insert capillary tubes into the correct wells of the heating block. 

(5) Bend the leads of the thermocouples under the remote cell гіп. 

(6) Insure that the zero temperature ice bath reference thermocouples 
are fully inserted into their glass tubes, which in turn are im- 
mersed in the ice and distilled water solution at 0°C. 

(7) Place a thermogram grid on the X-Y plotter. 

(8) Set the pen to room temperature on the T-scale. 

(9) Set the DTA 900 upper panel knobs: 

(a) T-zero shift. 

(b) T—scale, °C/in. 

(с) AT-zero shift. 

(d) AT-scale, °C/in. 

(e) Baseline slope. 

(10) Set the DTA 900 lower panel dials: 

(a) Temperature heating rate, range 0-30%C/min. 

(b) Set starting temperature dial to minimum setting to get the 

minimum voltage reading on the heater voltage indicator. 

(11) Set the program mode to HEAT. 

(12) Push RESET to erase previous program memory. 

(13) Place recorder pen to the DOWN position. 

(14) Place power knob to RECORD. 

(15) Start timer and Varian recorder. 

Figure 9 shows a sample thermogram for a heating rate of 20°C/min. 
beginning at room temperature. The endotherm at 204°C shows transfor- 
mation to the liquid phase, followed immediately by the exothermic de- 
composition with the peak in the exotherm at 251°C. Figure 10 shows a 
thermogram for ammonium nitrate. The endotherm at 40°C shows a transi- 
Elen from rhombic I to rhombic Il. At 91°C, a crystalline change to 


tetragonal occurs, then to cubic at 130°C. Melting occurs at 174°C 


followed by decomposition in the liquid phase. This run was made to 
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check the calibration of the machine, to give an indication of the 
sensitivity of the analyzer, and for use in a method of estimating the 
heat ОЁ reaction of RDX, to be described in the recommendations. The 
thermogram was compared with the ammonium nitrate thermogram presented 
on page 5-7 of the DTA Instruction Manual [Ref. 10]. The comparison 
showed a disparity in temperature for the peak exotherm, where that of 
Figure 10 showed a temperature approximately 10°C higher than that of 
the reference thermogram. This was corroborated by the Varian recorder, 


and subsequent temperature readings were taken from the Varian recorder. 


D. ISOTHERMAL EXPERIMENT 

A time base adaptor for the DIA 900 provides for maintaining the 
heating block at set temperatures over long periods of time. The X-Y 
plotter now provides for time in minutes on the abscissa and as before, 
AT on the ordinant. Experimental runs were made using the four mil- 
limeter test tube macro heating block. A variation of a method used 
in determining the effects of thermal aging of propellants was followed 
[Kuletz and Pakulak 1962]. This technique is known as microinmersion 
autocatalytic analysis (MAA). In this work, the sample of RDX was 
dropped into the well of the macro heating block after the block had 
been established at a desired temperature. The procedure followed is 
listed below. 

EXPERIMENTAL PROCEDURE FOR ISOTHERMAL EXPERIMENT 

(1) Turn DTA 900 power knob to STANDBY, allow one to two hours warm up. 


(2) Place inert reference material in both 4mm test tubes, tap and 
measure to a depth of four millimeters. 


(3) Place metal thermocouple holder caps over test tubes and insert 
thermocouples, insuring complete immersion. 


ΕἾ 





(4) Insert reference and sample 4mm test tubes in the large wells of 
the macro heating block. ‘Two small capillary tubes are used for 
the controlling and Varian recorder thermocouples. 

(5) Bend the thermocouple leads under the remote cell rim. 

(6) Check the zero temperature ice bath and thermocouples. 

(7) Place a thermogram grid on the X-Y plotter. 

(8) Set the time base adaptor to TIME. 

(9) Set the DTA 900 upper panel knobs: 

(a) T-zero shift. 

(о) T-scale, Min./in. 

(ο) AT-zero shift. 

(d) AT-scale. 

(e) Baseline slope. 
(10) Set the DTA 900 lower panel dials: 

(a) Set temperature heating rate at zero. 

(b) Set the starting temperature dial according to the scale of 

Figure 3G of section 3, DTA Instruction Manual [Ref. 10]. 

(11) Set the program mode to ISOTHERMAL. 


(12) Allow the temperature to stabilize at the desired setting and 
record the Varian reading. 


(13) Place a sample test tube of RDX in the macro heating block and 
reset the corresponding thermocouple. (Wear protective shield) 


(14) Push RESET to start zero time at the left end of the scale. 
(15) Place recorder pen to the DOWN position. 

(16) Record time and start timer. 

(17) Monitor Varian recorder. 


The resulting thermogram plots AT versus time. 
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V. RESULTS AND DISCUSSION 


The series of thermograms obtained from the variable heating rate 
method using RDX samples are reproduced in Figures 11-17. During each 
of these runs, temperatures were recorded at 5°C intervals and at sig- 
nifieant points on the thermogram (melting, peak exotherm) with the 
corresponding time in seconds. This data is plotted in Figures 198-36. 
An nt^ order regression analysis program of the WANG 700 was used to 
calculate the least squares fit to the data points, the slope of which 
is the actual heating rate for that run. The resulting equation gives 
the temperature at any time during each run for the programed heating 
rate. The calculated heating rates, along with the sample temperature 
of the peak in the exotherm Tm, are the required data for determining 
the kinetic parameters of the decomposition by Kissinger's method [Ref. 
2]. These data are summarized in Table I. 

The kinetic data for the decomposition reaction of RDX were deter- 
mined from the plot of the variable heating rate experimental data of 
Table I. The natural logarithm of the rate factor #/Tá was plotted 
versus the reciprocal of the Kelvin temperature of the heating block; 
as shown in Figure 37. The WANG 700 was used to calculate a least 
squares fit to the experimental data points. The slope of this line 
determines the activation energy which was calculated to be 49.4713 
keal/mole. The Arrhenius frequency factor was determined from Equation 
13, averaged over the temperature range of the experiment. The average 
value was 6.9199 x 1020, which varied from the fastest heating rate 
value by 1.1501 x 102° and from the lowest by 0.1415 x 1029. These 


values compare favorably with those of previous works [Rauch and 
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Wainright 1969, Robertson 1948}. The kinetic parameters are summarized 
in Table II. There is a 1.6% error in the value for the activation en- 
ergy when compared with that of Rauch and Wainright's work for Picatinny 
Arsenal. 

Appendix A is an IBM 360 computer readout which tabulates values for 
the heating rate with incremental increase in the peak temperature Tm. The 
experimentally determined values of activation energy and Arrhenius fre- 
quency factor were used in Equation 13 to generate values of the heating 
rate $. Appendix B tabulates the same information for comparison using 
the kinetic data determined by Rauch and Wainright. This data is presented 
in Table III along with values from the experimental runs for comparison. 
At higher heating rates, the experimental and calculated data differ toa 
larger extent. This is probably due to experimental time lag at the higher 
heating rates which caused lower temperature readings at the higher rates. 
Although the values of activation energies compare favorably, Table III 
shows that the experimentally determined temperatures of the peak in the 
exotherms at the various heating rates would have to be approximately 30°C 
higher to give heating rates comparable to the Picatinny data. 

The rate constant k over the experimental temperature range was calcu- 
lated using the Arrhenius equation. Appendix C tabulates values of k using 
the experimentally determined values of E4 and A. Appendix D presents the 
same data using the kinetic data determined in the study made for Picatinny 
Arsenal. These values are compared in Table IV. The smaller values for 
the rate constant resulting from the Pieatinny data correspond with the 


lower heating rate values shown in Table III. 
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Some variation existed in the peak temperature of the exotherms as 
evidenced in Table I. This resulted in the scatter noted in Figure 37. 
An attempt to reproduce the реак temperature at a given heating rate 
was made over three consecutive runs using the standard experimental 
technique. Results аге shown in Figure 38, where the three runs were 
carried out at a heating rate of 6%C/min. The maximum variation in the 
peaks was 4°C. This compares with the variation in temperatures for 
the experimental and calculated values of T, shown in Table ІІІ. The 
maximum deviation approaches .8% error, Senn is well within the accur- 
of the experiment. There is a definite increase in T, with heating 
rate as proposed by the DTA kinetic theory [Kissinger 1957]. The ef- 
feet of sample size appeared to be only an increase in the quantity of 
heat liberated, as shown in Figure 39. The peak temperature did not 
vary to any greater degree than in the normal manner with equivalent 
volumes. of material. 

The isothermal data gave somewhat less conclusive results. Тһе 
thermograms, Figures 40-44, show a difinite exothermic decomposition 
reaction which occurs at a constant block temperature after an induc- 
tion time which decreases with increased temperature. A temperature 
range of 192-210°C was investigated. The time elapsed to initiation 
of the exothermic reaction, or induction time, and the time elapsed to 
the peak in the exotherm, the time-to-deflagration, were observed and 
recorded. Table V summarizes the isothermal data. It was noted that 
although these reactions occur at or below the normal melting tempera- 
ture they occurred in the liquid phase. A physical observation during 
one of the runs confirmed a liquid phase decomposition with decomposi- 


tion gases bubbling off. The sample at 192°C showed no evidence of 
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reaction after six and one-half а Near the melting temperature 
(204°C) and above, the sample changed to the liquid phase immediately 
upon reaching block temperature with the exothermic decomposition en- 
suing. The microimmersion autocatalytic analysis (MAA) method was used 
to determine kinetic data from the isothermal phase of the experiment. 


This method is based on the equation; 


Т2. Аша - | 1 
- Se eEa/RT (15) 





where 09 15 the heat of reaction. A plot of the natural logarithm of 
the rate factor (1/T?) versus the reciprocal of the constant block tem- 
perature gives a slope equal to E,/R. In addition, an Arrhenius plot 
using the reciprocal of the time data versus 1000/T was made [Kuletz 
and Pakulak 1962]. Both induction time and time-to-deflagration data 
КЕ оса ае УІ lists Che iAuction time data. Two plots were 
made using these даса. Тһе natural logarithm of the reciprocal of the 
induction time versus the reciprocal of the isothermal block tempera- 
ture in degrees Kelvin is presented in Figure 45. The MAA procedure, 
plotting the natural logarithm of the rate factor (t./T?) versus 1000/ 
ТӨК, із shown in Figure 46. This was done again using the time-to- 
deflagration. The data for time-to-deflagration are listed in Table 
VII. The plots of these data are shown in Figures 47 and 48 respec- 
tively. The results are shown in Table VIII. In all cases the WANG 
{00 was used to determine a least squares fit to the data points and 
the resulting slopes. It can be seen that both the reciprocal time 
and MAA rate factor methods produce nearly identical results. The 


values obtained for E, are far too high for activation energies. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Experimental work done on HMX [Suryamarayana and Graybush 1967] 
suggested that materials which exhibit sigmoid pressure-time curves 
generally involve three stages in the decomposition reaction. Тһе 
first stage is the induction period. This is followed by an acceler- 
ation period where the rate reaches a maximum. The material is con- 
sumed in the decay period. Even with this mass spectrometric study, 
the conclusions were limited due to the sagan decomposition products 
formed. It was observed that each product reached a constant rate of 
formation after an acceleratory stage which appeared due to autocat- 
alysis. It was found that HMX exhibited different modes of decomposi- 
tion. The kinetic data changed in different discrete temperature 
ranges. These findings were supported in a DTA analysis of HMX [Hondee 
1971], where three separate activation energies were found during three 
discrete stages of the decomposition process. 

Causes for acceleration of the decomposition of HMX to a constant 
rate can be attributed to [Suryamarayana and Graybush 1967]: 

(1) Progressive melting as a result of lowered melting point by 

mime decomposition products. 

(2) Self-heating. 

(3) Autocatalysis by products. 

(4) Acceleration due to structural factors such as the increase 

ра спе отер of nuclei. 

During this experimental study of the decomposition of RDX there 
was no evidence of different activation energies since the heating rate 
in all runs remained constant over the entire decomposition. The peak 
exotherm, although somewhat variable, did not exhibit more than one 
sharp peak. The thermograms for the variable heating rate experiments 


> 
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(Figures 11-17) show that the exothermic decomposition reached a peak 
which was in most cases followed by a second peak. This second peak 
15 probably a result of the capillary action of the liquid sample 
causae che thermocouple to rise in the tube, then drop back down just 
before the reaction completes causing another jump in the thermogram. 

The isothermal experiment thermograms (Figures 40-44) show a def- 
inite decomposition reaction after an induction time at lower tempera- 
tures. This is in agreement with previous. studies which have found 
evidence of self-heating and autoeatalysis in the decomposition of RDX 
and other explosives. These data can be applied to estimates of the 
sensitivity of explosives and propellants. 

Several reasons may be proposed for the discontinuities in the 
variable heating rate experiment. Although the RDX sample was assumed 
homogeneous and pure, it may well have been nonuniform in nature and 
composition. Possible water content in the samples was neglected. 

DTA results vary between laboratories due to the difference in experi- 
mental techniques and equipment as well as in the interpretation of 
the results. Standardization of method would be beneficial. 

Explosive decomposition rates vary rapidly with temperature and 
for ideal results the temperature of the sample should at all times 
bus. This is not the true situation due to thermal conduction 
in the sample and self-cooling or self-heating effects arising from 
the reaction. Robertson [Ref. 9] showed that errors in the values of 
the kinetic parameters occur due partly to autocatalysis and self- 
heating. Autocatalysis increases the activation energy value by ap- 


proximately 2000 cal/mole as a maximum error and causes a corresponding 
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variation in the Arrhenius frequency factor by as much as a factor of 
10. Errors due to self-heating are also incurred and tend to increase 
Es and Á as well. This effect increases with increased temperature. 

Experimental errors are also induced by the apparatus. The accu- 
racy of the DIA method depends on the precision with which a uniform 
heating rate can be maintained. Although the DIA 900 provides for a 
constant rate to the heating block, the exterior of the sample is hot- 
ter than the interior and reacts sooner and faster. This causes addi- 
tional autocatalytic and self-heating effects which are accented by an 
increased heating rate. Another source of error in this work was the 
uncontrolled atmosphere. This provides for possible variations in the 
data from day to day. Bubbles of the decomposition gases which formed 
en the thermaeouples caused the high intensity chatter on the thermo- 
grams. Capillary action of the melt during decomposition caused it 
to climb the walls of the tube surrounding the thermocouple sleeve. 
This drew the material away from the thermocouple probe. Care had to 
be taken to keep the level of the liquid below the sleeve in order to 
ΠΣΕ this effect. 

Notwithstanding these adverse and random effects, the dominant fac- 
tor in controlling the shape and position of the endotherms and exo- 
therms is the nature of the reaction itself. The favorable results 
obtained lend credibility to the DI'A method and suggest that further 
development in the technique will prove useful in providing a fast and 
reasonably accurate method of characterizing the cessent and kinetic 
properties of new explosives and propellants. 


Areas for further study are wide and varied.  Considerably more 


information as to the mechanisms of' the reaction can be obtained by 
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combining the DTA method with one or more of several more precise an- 
alytical techniques such as thermogravimetry, infrared spectrometry, 
mass spectrometry, X-ray measurements, and visual observations. 

More information may be determined from the DTA thermograms them- 
selves. It has been shown that the area under (һе peak in the thermo- 
gram is proportional to the heat of reaction [Rivette and Besser 1961]. 
This provides a means to quickly estimate the heats of reaction of the 
sample. The shape of the exotherm provides a means to estimate the 
reaction order [Kissinger 1957]. Rough estimations using this method 
showed the reaction order to vary from 0.2 to 0.8 over the range of 
heating rates considered. Further analysis of the isothermal thermo- 
grams appears warranted to discover the significance of the two slopes 
on the low temperature side of' the exotherm and the shape and extent 
of the decay pattern on the high temperature side, as well as the sig- 
nificance of the high values obtained for E, from the Arrhenius and 


MAA plots of the induction time and time-to-deflagration data. 
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ТАВГЕ 1 


EXPERIMENTAL DATA, VARIABLE HEATING RATE THERMAL ANALYSIS 


RUN # %,90/МІМ Tm °K in(4/%%) — Q/Tmyao3 
1 2.431 504.00 -11.5601 1.9802 
2 2.880 3000 -11.3754 1.9960 

3.333 509.00 -11.2610 1.9646 

4 3.8571 ο ο -11.1306 1.9493 
5 4.3636 | 505525 -10.9768 1.9792 
6 4.6957 509.00 -10.9183 1.9646 
T 5.1940 520250 -10.8233 1.9589 
8 5.6552 514.50 -10.7538 1.9436 
9 6.2069 515.25 -10.6636 1.9408 
10 6.8182 513.00 -10.5610 1.9493 
11 7.9149 51.25 -10.4283 1.9333 
12 8.8696 POTS -10.2891 1.9579 
13 9.7143 511,50 -10.2011 1.9550 
14 10.5556 512,75 -10.1229 1.9503 
15 11.6471 519.00 ~10. 0488 1.9268 
16 16.5217 517200 -9.6914 1.9342 
17 21.0000 524.50 -9.4804 1.9066 
18 26.0769 55:25 -9.2591 1.9111 
19 30. 0000 528.00 -9.1370 1.8939 
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TABLE 11 


KINETIC PARAMETERS 


38 


Ea (kcal/mole) А (ѕес-1) 
EXPERIMENTAL Ң9.5 επ, ΙΡ!» 
RAUCH & WAINRIGHT 48.7 ПОО 
ROBERTSON 47.5 er sp 5 
TABLE III 
HEATING RATE DATA COMPARISON 
PEAK 
TEMPERATURE HEATING RATE ф 
_ то EXPERIMENTAL CALCULATED PICATINNY 
504.00 2.4310 2.4817 0.1338 
509.00 4.6957 4.1122 0.2200 
513.00 6.8182 6.1167 0.3253 
524.50 21.0000 18.5317 0.9694 
528.00 30.0000 25.7247 1.3391 
TABLE IV 
RATE CONSTANT DATA COMPARISON 
PEAK 
TEMPERATURE RATE CONSTANT k (sec”“) 
_ то EXPERIMENTAL PICATINNY 
504.00 0.2433 0.0129 
509.00 0.3952 0.0208 
513.00 0.5787 0.0303 
524,50 1.6772 0.0864 
528.00 2.2974 ὉΠ. 








> 


о ἈΠ 


TABLE V 


EXPERIMENTAL DATA, ISOTHERMAL ANALYSIS 


ATmax 


· 39 
250 
.66 
1250 
1.50 
1.60 


3.20 


| INDUCTION TIME DEFLAGRATION 
T (CC) T CK) to (min) to (sec) т (піп) q (sec) 
192.00 165.00 EXCESS OF 6-1/2 HRS. 
MOUS Шбфесо 092.0 12,520 9238.5 14310 
196.75 469.75 130.0 7,848 145.6 8,736 
ПОО 471.00 85.0 5.1000, ΠΡ 6,450 
198.50 171.50 37.5 2,250 58.5 3,510 
199.50 472.50 15.0 900 34.8 2,088 
204.50 177.50 10.8 648 19.2 1,152 
206.75 179.75 0.0 0 0.0 0 
TABLE VI 
ISOTHERMAL ANALYSIS, INDUCTION TIME DATA 
1000/f to sec 1/to ας E E) 
2.1402 11,520 .8681x10* -9.3518 -2.9419 
2,1288 7,8488 1274103 -8.9680 -3.3364 
21251 5,100 .1961х103 -8.5370 ΠΠ 
2.1209 22500 ЩОБ < т ТА -4.5932 
2.1161 9000 .1111х102 -6.8024 - 
2.0964 648 — .1503x10 ^ -6.4739 
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ISOTHERMAL ANALYSIS, TIME-TO-DEFLAGRATION 


1000/7 
2.1402 
2.1288 
2.1231 
2.1209 
2.1164 


2.0964 


TABLE VII 


1/τ(5δες” 1) 
.6988x107* 
1145x107? 
1550x107? 
. 2819x107? 
. 4789x107? 
.8681x107? 


TABLE VIII 


ISOTHERMAL RESULIS 


INDUCTION TIME 


ios m 


-6.8194 
(290239 
-7:6322 


-ὂ. 2128 


-8.7665 
-9.3801 


1n(1/1) 
-5.4744 
-4.9809 
-4.6775 
-4.0690 
-3.5496 
-2.9549 


TIME-TO-DEFLAGRATION 


Ea (kcal/mole)> 159.0282 160.9012 


ДО 


131.0144 132.8778 
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Figure ІІ. Variable Heating runs 1-3. 
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0.296 


heating rate = 4.36% C/min. 
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Figure l2. Variable Heating rate runs 4&5. 
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heating rate = 4./0%C/min. 


heating rate = 5.19% C/min. 


heating rate = 5.66? C/min. 


| 


204 220 240 260 280 300°C 
Figure I3. Variable Heating Rate Runs 6-8. 
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heating rate = 6.21% C/min. 


0.5° C 


heating rate = 6.82° C/min. 
ος 


heating rate = 1.92% C/min. 


| 


204 220 240 260 280 οσο ఈ 
Figure I4. Variable Heating Rate Runs 9- Il. 
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heating rate = 8.87 ° C/min. 


heating rate = 9. 71° C/min. 


heating rate =0.56* C/min. 


ROC 


| 











Πλ “προ οσο 200 280 300°C 
Figure 15. Variable Heating Rate Runs 12-14. 
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heating rate 211.65? C/min. 


heating rate 216. 52? C/min 





heating rate 221.00? C/min 


КӨС 


204 220) 240 260 280 TITAN 
Figure l6. Variable Heating Rate Runs 19-17. 
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heating rate = - 
26.О8° С/т\їп. 


heating rate= 
3O.OO° C/min. 


LOSE 


204 220 240 260 280 SOON 
Figure If. Variable Heating Rate Runs 18 8 19. 
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APPENDIX А 


HEATING RATE VS. PEAK TEMPERATURE, 500.00-530.509K 
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APPENDIX B 


HEATING RATE VS. PEAK TEMPERATURE, 500.00-530.509К 
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APPENDIX C 


RATE CONSTANT VS. ISOTHERMAL TEMPERATURE, 470.00-550.00 9К 
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APPENDIX C (continued) 
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APPENDIX D 


RATE CONSTANT VS. ISOTHERMAL TEMPERATURE, 470.00-550.00 °K 
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APPENDIX D (continued) 
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